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The  flow-field  for  reactant  distribution  is  an  important  design  factor  that  influences  the  performance 
of  polymer  electrolyte  membrane  fuel  cells  (PEMFCs).  Under-rib  convection  between  neighboring  chan¬ 
nels  has  been  recognized  to  enhance  the  performance  of  PEMFCs  with  serpentine  flow-fields.  This  study 
presents  a  simple  design  method  to  generate  multi-pass  serpentine  flow-fields  (MPSFFs)  that  can  maxi¬ 
mize  under-rib  convection  in  a  given  cell  area.  Geometrical  characterization  indicates  that  MPSFFs  lead 
to  significantly  higher  under-rib  convection  intensities  and  more  uniform  conditions,  such  as  reactant 
concentrations,  temperature,  and  liquid  water  saturation,  compared  with  conventional  serpentine  flow- 
fields.  The  implications  of  the  enhanced  under-rib  convection  due  to  MPSFFs  on  the  performance  of 
PEMFCs  are  discussed. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  produce  fewer  harmful  emissions  and  are  more  effi¬ 
cient  compared  with  heat  engines  (Carnot  efficiency)  by  directly 
converting  the  chemical  energy  of  fuels  into  electricity  without 
combustion  [  1,2 ].  In  addition,  fuel  cells  can  overcome  the  drawbacks 
of  present  rechargeable  batteries  by  providing  longer  operation 
time  and  instant  refuelling  features.  Thus,  fuel  cells  are  expected 
to  serve  as  important  power  sources  in  many  applications  in  the 
near  future.  Among  the  many  types  of  fuel  cell,  polymer  electrolyte 
membrane  fuel  cells  (PEMFCs)  have  been  applied  to  a  diverse  vari¬ 
ety  of  uses  that,  include  portable,  automobile  and  stationary  power 
generation.  The  relatively  low  operating  temperature  of  PEMFCs, 
i.e.  below  100  °C  enables  fast  start-up  at  the  beginning  phase  and 
good  dynamic  response  to  load  change.  These  excellent  transient 
characteristics  of  PEMFCs  are  essential  for  applications  that  require 
short-term  repeated  operations.  Currently,  much  effort  is  being 
applied  to  develop  more  efficient  PEMFCs  that  are  able  to  operate 
at  higher  current  densities. 

The  performance  of  a  PEMFC  is  primarily  determined  by  the 
intrinsic  electrochemical  efficiency  of  the  membrane  electrode 
assembly  (MEA).  Nevertheless,  other  factors,  such  as  flow-field 
design,  thermal  and  water  management,  and  operational  control 
are  also  important  [3].  Accordingly,  many  workers  have  investi¬ 
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gated  viable  flow-field  designs  suitable  for  PEMFCs  [4].  Three  basic 
type  of  flow-field  are  commonly  used  for  PEMFCs,  namely,  paral¬ 
lel,  serpentine  and  interdigitated  flow-fields.  The  parallel  version 
is  the  simplest  and  requires  the  smallest  pressure  drop  by  equally 
distributing  the  flow  rate  into  many  straight  parallel  paths.  On  the 
other  hand,  non-uniform  distribution  of  the  reactant  gas  between 
paths  can  occur  in  parallel  flow-fields  if  the  flow  resistance  of  each 
path  is  not  maintained  at  a  comparable  level.  Serpentine  flow-field 
designs  have  been  proposed  to  solve  the  maldistribution  problem 
inherent  to  parallel  flow-fields  by  forcing  the  reactant  gas  to  flow 
through  a  single,  long,  meandering  path  that  travels  over  the  entire 
active  area.  Thus,  the  flow  speed  and  the  pressure  drop  in  serpentine 
flow-fields  are  both  very  high,  which  also  facilitates  a  two-phase 
transport  of  liquid  water  droplets  through  the  gas  channels  of  PEM¬ 
FCs. 

Interdigitated  flow-fields  are  based  on  dead-end  channel 
designs,  which  direct  reactant  gases  to  flow  through  porous  gas 
diffusion  layers  (GDLs)  by  under-rib  paths  [5,6].  The  presence  of 
a  convective  flow  in  the  under-rib  regions  enables  more  effective 
utilization  of  electrocatalysts  by  increasing  reactant  concentration 
and  facilitating  liquid  water  removal  in  those  regions.  This  under¬ 
rib  convection  has  recently  been  recognized  by  many  researchers 
[7-17]  as  a  non-negligible  transport  process  that  influences  the 
performance  of  PEMFCs  and  direct  methanol  fuel  cells  (DMFCs) 
with  serpentine  flow-fields.  Experimental  studies  have  shown  that 
higher  GDL  permeability  improves  the  performance  of  PEMFCs  with 
serpentine  flow-fields  [8,9]  and  reduces  the  pressure  drop  [16]. 
Pharoah  [10]  observed  that  under-rib  convection  can  not  be  ignored 
when  GDL  permeability  exceeds  10-13  m2,  which  is  consistent  with 
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the  analytical  model  prediction  for  the  relative  influence  of  under¬ 
rib  convection  [13].  The  role  of  under-rib  convection  in  DMFCs  was 
also  numerically  studied  by  Ye  et  al.  [14]. 

Recently,  Xu  and  Zhao  [18]  proposed  a  new  flow-field  design 
for  polymer-based  fuel  cells,  namely,  the  convection  enhanced  ser¬ 
pentine  flow-field  (CESFF).  They  experimentally  showed  that  a 
small  DMFC  (3  cm  x  3  cm)  with  a  CESFF  could  perform  better  and 
more  stably  than  a  similar  DMFC  with  a  conventional  serpentine 
flow-field,  especially  at  smaller  air  flow  rates.  They  attributed  the 
enhanced  performance  to  the  ability  of  the  CESFF  to  handle  liquid 
water  in  the  under-rib  regions.  Similar  flow-field  designs  have  also 
been  proposed  as  coolant  flow-fields  for  PEMFCs  [19,20].  By  plac¬ 
ing  flow  channels  of  different  path-lengths  (the  travelled  distance 
from  an  inlet)  near  one  another,  an  improved  temperature  unifor¬ 
mity  in  the  cooling  plate  was  numerically  predicted.  It  should  be 
noted  that  uniform  temperature  is  also  favourable  for  the  operation 
control  and  the  long-term  durability  of  PEMFCs.  The  idea  of  close 
placement  of  flow-channels  with  different  path-lengths  was  previ¬ 
ously  proposed  by  Kaufman  and  Terry  [21  ].  Later,  Qi  and  Kaufman 
[22]  demonstrated  that  a  PEMFC  with  a  double-path  (counter¬ 
flow)  cathode  flow-field  could  operate  stably  under  low-humidity 
conditions  when  using  dry  reactant  gases.  They  reported  that  the 
double-path-type  flow-field  facilitated  the  humidification  of  dry 
entering  gases  with  internally  produced  water. 

Based  on  previous  studies  [18-22],  multi-pass  serpentine  flow- 
fields  (MPSFFs)  are  proposed  in  this  study  to  enhance  the  under-rib 
convection  and  thus  improve  the  performance  of  PEMFCs.  The 
MPSFFs  cause  reactant  gases  to  travel  mutiple  timer  between  the 
region  near  an  inlet  and  that  near  an  outlet  and  thereby  results  in 
closely-interlaced  channels  with  widely  different  path-lengths.  The 
channel  patterns  in  MPSFFs  are  believed  to  enable  higher  under¬ 
rib  convection,  more  uniform  concentration  distribution,  and  fast 
exhaust  of  liquid  water  from  the  under-rib  regions.  This  study 
presents  a  step-by-step  design  to  generate  MPSFFs,  along  with 
several  exemplary  flow-fields,  i.e.,  single-path  MPSFFs  for  a  small 
active  area  of  9  cm2  (3  cm  x  3  cm)  and  parallel-path  MPSFFs  for  a 
relatively  large  active  area  of  81  cm2  (9  cm  x  9  cm).  Next,  the  geo¬ 
metrical  characterization  is  presented  for  those  exemplary  single- 
and  parallel-path  MPSFFs  and  is  based  on  the  channel  path-length 
distribution  in  the  active  area.  Finally,  the  effects  of  under-rib  con¬ 
vection  enhancement  by  MPSFFs  on  the  overall  performance  of 
PEMFCs  are  discussed  in  relation  to  influence  of  cell  size. 


2.  Theory  and  calculations 

2.1.  Under-rib  convection  mechanism 


The  mechanism  of  under-rib  convection  is  explained  with  three 
neighbouring  channel  regions  A,  B,  and  C  in  the  exemplary  (3-pass) 
serpentine  flow-field  shown  in  Fig.  1(a).  The  path-lengths  ZA,  ZB, 
and  Zc  are  defined  as  the  distance  travelled  by  the  reactant  gas 
from  an  inlet  to  the  three  regions  along  the  channel— illustrated 
by  dotted  lines  in  Fig.  1(a).  The  path-length,  z,  can  be  viewed  as  a 
coordinate  whose  upper  limit  is  the  same  as  the  total  length  of  the 
path,  zmax.  The  pressure  and  the  reactant  concentration  in  PEMFCs 
are  expected  to  change  linearly  along  the  flow  channel,  i.e.,  with 
respect  to  z. 

The  total  flow  rate,  qtot,  depends  on  the  total  pressure  drop, 
Aptot.  across  the  flow-field  as  espressed  by  [23]: 


Qtot  — 


Kch^  Agtot 

P'  Zmax 


(1) 


where:  /x  is  the  viscosity;  Ach  is  the  cross-sectional  area  of  the 
channel  (height  hch  x  width  wch);  Kch  is  the  channel  permeabil¬ 
ity  derived  from  the  Flagen-Poiseuille  equation.  Note  that  Eq.  (1) 
assumes  that  the  flow  rate  is  constant  along  the  channel  with  a  low 
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Fig.  1.  Under-rib  convection  and  transport  in  serpentine  flow-field:  (a)  definition  of 
path-length,  z;  (b)  under-rib  convection  and  transport  mechanisms. 


under-rib  convection  flow  rate  (qrib^0)»  and  that  the  secondary 
(non-linear)  pressure  loss  in  90°  or  180°  corners  is  negligible.  Thus, 
the  permeability,  Kch,  is  determined  by  [24]: 


K ch  = 


Jk. 

28.45’ 


(2) 


from  the  Darcy  friction  factor  of /=  56.9/Re  for  square  channels.  In 
this  study,  Kch  is  estimated  to  be  3.515  x  10-8  m2  for  square-channel 
geometry  with  hch  =  wch  =  1  mm.  The  pressure  difference  between 
two  neighbouring  channel  regions  is  expressed  in  terms  of  the  path- 
length  difference  as: 


IPa-PbI 


Aptot 

Zmax 


\ZA  -  ZB I 


gtotfL 

^ctAh 


\za  -zb|  • 


(3) 


Let  us  assume  that  the  flow-field  shown  in  Fig.  1(a)  is  for  the 
cathode  of  an  PEMFC.  The  reactant  concentration,  c  (mol  m-3 ),  then 
denotes  the  molar  concentration  of  oxygen,  and  the  total  concen¬ 
tration  loss,  Actot  =  Qn  -  Cout,  across  the  flow-field  is  calculated  from 
the  mass  balance  of  oxygen  as  according  to: 


(Cjn  —  C0ut)<?tot  —  ACtot  Qtot  — 


(4) 


Here,  /tot  (A)  is  the  current  generated  over  the  entire  cell  area  and 
F  is  the  Faraday  constant  (96487  C  mol-1 ).  Note  that  Eq.  (4)  ignores 
the  change  in  the  total  flow  rate,  qtot,  due  to  oxygen  consumption. 
In  fact,  the  reactant  flow  rate  in  the  cathode  is  relatively  constant 
because  the  humidified  air  (with  an  oxygen  mole  fraction  of  about 
0.15)  is  supplied  at  high  excess  air  ratios.  If  current  is  generated 
uniformly  throughout  the  active  area,  the  concentration  difference 
between  two  neighbouring  channel  regions  is  also  expressed  in 
terms  of  the  path-length  difference,  i.e., 


|cA-cB| 


Actot 

Zmax 


|zA-zb| 


hot 

4FqtotZmax 


\Za  -  ZB I 


(5) 


Fig.  1(b)  shows  the  cross-sectional  view  of  the  three  adjacent 
channel  regions  A,  B,  and  C.  Region  A  has  the  highest  pressure  and 
reactant  concentration,  while  region  C  has  the  lowest  among  the 
regions  (pA  >pB  >Pc  and  cA  >  cB  >  cc).  Due  to  the  pressure  difference 
between  neighbouring  regions,  convective  flow  occurs  through 
under-rib  regions  in  the  GDL,  as  shown  in  Fig.  1(b).  This  additional 
convective  flow  through  GDLs  is  called  the  under-rib  convection, 
and  is  believed  to  improve  the  performance  of  PEMFCs.  Under-rib 
convection  increases  the  reactant  concentration  in  the  under-rib 
regions,  facilitates  liquid  water  removal  from  those  regions,  and 
enables  a  more  uniform  concentration  distribution  throughout  the 
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flow-fields.  The  under-rib  convection  flow  rate,  grib,  is  proportional 
to  the  pressure  difference  between  two  neighbouring  regions.  For 
example,  qA®  between  regions  A  and  B  in  Fig.  1(b)  can  be  expressed 
as  [23]: 


„AB  _  ^gdl  aab  IPa  -  Pb  I  _  „  ^gdl^rib  \Z&  -  Zg  | 
q"h  -  H  rib  wrib  -qtotI<chAch  wrib 


(6) 


where  Kgdl  is  the  flow  permeability  of  the  GDL,  wrib  is  the  rib  width, 
and  A^®  is  the  flow  area  for  the  under-rib  convection.  Note  that 
AA®  is  determined  by  multiplying  the  GDL  thickness,  tgdl,  and  some 
length,  Zch,  of  the  regions  along  the  channel  (the  under-rib  convec¬ 
tion  through  the  catalyst  layer  is  ignored). 

Fig.  1(b)  shows  that  other  transport  mechanisms  also  become 
active  when  channel  regions  that  have  different  path-lengths  are 
in  close  contact.  The  local  temperature  and  liquid  water  satura¬ 
tion  in  GDLs  depend  on  the  electrochemical  reaction  rate,  and  this 
is  generally  governed  by  the  reactant  concentration,  which  varies 
along  the  channel  (with  respect  to  z).  The  temperature  gradient  due 
to  closely  placed  channels  of  different  path-lengths  facilitates  con¬ 
ductive  heat  transfer  by  the  under-rib  or  cross-rib  paths,  and  the 
water  saturation  gradient  facilitates  capillary  transport  of  water  by 
the  under-rib  path.  These  transport  mechanisms  may  be  termed  as 
under-rib  transports  to  distinguish  them  from  under-rib  convec¬ 
tion.  While  under-rib  convection  occurs  in  only  one  direction  (from 
a  high-pressure  region  towards  a  low-pressure  region),  under-rib 
transports  may  occur  in  both  directions,  as  shown  in  Fig.  1(b).  In 
fact,  Qi  and  Kaufman  [22]  attributed  the  successful  dry  operation 
of  their  PEMFC  to  the  under-rib  transport  of  water  in  the  opposite 
direction  of  the  under-rib  convection.  Thus,  under-rib  convection 
and  transport  are  believed  to  enhance  the  uniformity  of  operation 
conditions  in  the  flow-fields,  including  the  reactant  and  product 
concentrations,  temperature  and  liquid  water  saturation,  as  shown 
in  Fig.  1(b). 


2.2.  Geometrical  characterization 

The  two  unit  cells  used  for  the  generation  and  characterization 
of  flow-fields  are  presented  in  Fig.  2(a).  The  unit  cell  length,  Zunit, 
is  defined  as  the  sum  of  the  channel  width,  wch,  and  the  rib  width, 


wrib,  in  Fig.  1  (b).  Thus,  Zunit  is  2  mm  in  this  study,  while  both  wch  and 
wrib  are  assumed  to  be  1  mm.  It  should  be  noted  that  any  serpentine 
flow-field  for  PEMFCs  can  be  constructed  by  paving  the  active  area 
with  channel  cells  that  have  been  chosen  from  the  two  unit  cells 
and  then  rotated  properly. 

A  channel  cell  generally  has  two  neighbouring  channel  cells  con¬ 
nected  through  the  under-rib  regions,  which  is  indicated  by  the 
double-headed  arrows  in  Fig.  2(a).  The  channel  cells  located  at  the 
edge  of  the  flow-field  may  have  no  or  only  one  neighboring  channel 
cell.  The  path-length  difference,  Azit  of  a  unit  cell  i  is  determined 
by  summing  the  absolute  difference  of  the  path-length  z  between 
the  cell  i  and  each  neighbouring  cell  j,  or 


AZi  =  \  Z[  —  Zj  (7) 

z — ^neighboring  channels  ofi  1  1 

The  calculation  procedures  for  A  zx  of  the  channel  cells  A,  B,  and  C 
are  shown  in  Fig.  2(b).  The  Azj  is  a  good  geometrical  index  for  mea¬ 
surement  of  the  under-rib  convection  in  a  flow-field,  as  explained 
by  Eq.  (6). 

The  local-averaged  path-length,  z,  described  in  Fig.  2(b)  is 
another  good  geometrical  index  for  measuring  the  uniformity  of 
local  conditions  in  a  flow-field.  Conditions  such  as  the  reactant  and 
product  concentrations,  temperature  and  liquid  water  saturation, 
generally  vary  along  the  channel  (with  respect  to  z).  The  zx  of  a  unit 
channel  cell  i  is  calculated  by  averaging  its  path-length  with  those 
of  adjacent  channel  cells  j.  A  3  x  3  averaging  mask  (a  dotted  box  in 
Fig.  2(b))  consisting  of  nine  adjacent  cells  is  used  to  determine  zx 
as: 


—E  zj 

navg  Z- — q'=adjacent  channels  ofi 


(8) 


where  naVg  is  the  number  of  adjacent  cells  used  for  the  local  aver¬ 
aging.  Thus,  nav g  is  nine  for  most  channel  cells,  but  for  channel  cells 
located  at  the  edge  of  the  flow-field,  a  navg  of  six  or  four  should  be 
used,  as  shown  in  Fig.  2(b). 

A  small  standard  deviation,  oz ,  of  the  local-averaged  path- 
length,  z,  from  its  average  of  (z)avg  =  zmax/2,  is  desirable  to  obtain 
more  uniform  operational  conditions  in  the  entire  active  area.  In 
addition,  the  distribution  of  z  can  be  used  to  estimate  the  perfor¬ 
mance  of  cooling  plates  for  PEMFCs  because  the  temperature  of  a 
cooling  fluid  increases  rather  linearly  with  z  [19,20]. 


ZA  =(z1  +  z2  +  23  +  z5  +  zA+z6  +  Z7  +  ^  +  29)/9 

ZB  =(%+z4  +  z6+zB+z9  +  zc)/6 
ZC  =(Z6+ZB+Z9  +  Zc)/4 


Fig.  2.  Geometrical  characterization  of  flow-fields:  (a)  unit  channel  cells;  (b)  defi¬ 
nitions  of  path-length  difference,  A z,  and  local-averaged  path-length,  z. 


2.3.  Under-rib  convection  intensity 

The  under-rib  convection  intensity  can  be  assessed  by  com¬ 
paring  the  under-rib  convection  flow  rate,  qrib,  and  the  total  flow 
rate  across  the  flow-field,  qtot.  This  approach  has  been  previously 
used  to  evaluate  the  under-rib  convection  in  serpentine  flow-fields 
[10-13,18].  In  the  MPSFFs  proposed  in  this  study,  the  first  one-third 
of  the  path-length  of  the  channel  region  (z<zmax/3)  generally  con¬ 
tacts  with  a  channel  region  with  a  larger  path-length  (z>zmax/3). 
Thus,  the  reactant  flow  rate  along  the  channel  qch  almost  always 
decreases  forz<zmax/3  due  to  the  under-rib  convection  in  MPSFFs. 

In  this  study,  the  under-rib  convection  intensity  is  defined  as 
the  ratio  of  the  under-rib  convection  flow  rate,  qribl/3,  for  the  first 
one-third  of  path-length  (z<zmax/3)  to  the  total  flow  rate,  qt0 t-  By 
setting  Arib  in  Eq.  (6)  as  zmaxtgdl/3,  the  under-rib  convection  inten¬ 
sity  Qrib.i /3 / Qtot  is  determined  as: 

to,  1/3  _  1  ^gdl£gdlzmax  (^)avg,l/3 

qtot  ”  3  JCchhchwch  wrib  ’ 

where  (A z)avgl/3  is  the  average  path-length  difference  for  the  first 
1/3  path-length  of  the  flow-field.  Eq.  (9)  clearly  shows  that  the 
under-rib  convection  intensity  depends  on  the  maximum  path- 
length,  zmax,  and  therefore  indicates  that  cell  size  is  an  important 
factor  in  determining  the  under-rib  convection  intensity.  Although 
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Table  1 

Geometrical  parameters  for  generating  MPSFFs  and  estimating  their  under-rib  con¬ 
vection  intensities. 


Parameter 

Explanation 

Value 

hch 

Channel  height 

1  mm 

Wch 

Channel  width 

1  mm 

wrib 

Rib  width 

1  mm 

fgdl 

GDL  thickness 

0.2  mm 

Kch 

Channel  permeability 

3.515  x  10“8  m2 

^gdl 

GDL  permeability 

1.0  x  10“12  m2 

Eq.  (9)  allows  an  easy  comparison  of  the  under-rib  convection 
intensities  in  different  flow-fields,  it  only  gives  approximate  values 
because  a  small  under-rib  convection  is  assumed  in  its  derivation. 

Geometrical  parameters  for  estimating  the  under-rib  convec¬ 
tion  intensities  of  MPSFFs  are  summarized  in  Table  1.  The  GDL 
permeability,  JCgdl,  is  assumed  to  be  1  xl0_12m2  in  this  study. 
The  permeabilities  of  different  GDL  materials  for  PEMFCs  have 
been  measured  by  several  researchers  [8,25,26].  Williams  et  al. 
[8]  reported  the  through-plane  permeability  of  several  GDL  mate¬ 
rials,  and  obtained  8.7  x  10-12  m2  for  a  bare  Toray  carbon  paper 
(TGP-H-120)  with  75.6%  porosity,  and  1.9  x  10-12  m2  for  the  same 
paper  after  a  microporous  layer  application.  Similarly,  Feser  et  al. 
[26]  showed  that  the  in-plane  permeability  of  a  Toray  carbon  paper 
(TGP-H-60)  amounts  to  about  5  x  10-12  m2  for  70%  porosity.  Con¬ 
sidering  the  reduction  of  pore  volume  in  GDL  materials  due  to 
hydrophobic  coating  and  compression,  the  assumed  GDL  perme¬ 
ability  of  1  x  10-12  m2  seems  to  be  reasonable.  In  addition,  the 
presence  of  liquid  water  in  GDLs  during  PEMFC  operation  further 
reduces  the  gas-phase  permeability  due  to  two-phase  flow  effects. 

3.  Results 

3.1  Basic  channel  patterns 

Fig.  3(a)  depicts  several  basic  channel  patterns  denoted  accord¬ 
ing  to  the  number  of  passes,  npass.  AH  of  the  channel  patterns  except 

(a) 

1  -pass 

2 -pass 

3-pass  ■ 

I  1  5-pass  (spiral) 


4-pass 

1 


7-pass  (spiral) 

5-pass 

■  r"~ — — i 

i  ■ 

i  i _ 


Fig-  3.  Construction  of  basic  parallel  flow-fields:  (a)  basic  multi-pass  channel  pat¬ 
terns;  (b)  parallel  arrangement  of  3-pass  channel  patterns. 


for  the  1-pass  pattern  force  the  reactant  gas  to  travel  between 
the  inlet  region  and  the  outlet  region  multiple  times.  This  multi¬ 
pass  flow  structure  incurs  closely-interlaced  channels  of  widely 
different  path-lengths— which  are  essential  for  constructing  the 
multi-pass  serpentine  flow-fields  (MPSFFs)  proposed  in  this  study. 
Therefore,  the  basic  channel  patterns  shown  in  Fig.  3(a)  serve  as  the 
elementary  building  blocks  for  generating  single-  or  parallel-path 
MPSFFs.  Fig.  3(b)  shows  an  exemplary  parallel  flow-field  that  is  con¬ 
structed  by  arranging  the  3-pass  patterns.  Such  parallel  flow-fields 
are  considered  first  because  they  have  relatively  simple  geometries 
suitable  for  characterizing  the  basic  channel  patterns. 

The  1-pass  pattern  does  not  produce  a  multi-pass  flow  struc¬ 
ture  when  arranged  in  parallel,  but  instead,  serves  as  a  reference 
corresponding  to  conventional  parallel  flow-fields  with  negligible 
under-rib  convection.  The  2-pass  pattern  is  the  simplest  among  the 
multi-pass  channel  patterns,  and  is  also  expected  to  lead  to  the 
highest  under-rib  convection  intensity.  The  2-pass  pattern  does  not, 
however  easily  generate  parallel  flow-fields  because  the  inlet  and 
outlet  ports  must  be  made  on  the  same  side  of  the  cell  area.  In  fact, 
any  basic  channel  pattern  having  an  even  number  of  passes,  i.e.,  the 
2-  and  4-pass  patterns  in  Fig.  3(a),  is  not  suitable  for  constructing 
parallel-path  MPSFFs.  The  3-pass  pattern  is  the  simplest  among  the 
basic  patterns  with  odd  numbers  of  passes,  and  is  also  suitable  for 
both  single-  and  parallel-path  MPSFF  construction. 

The  4-pass  and  5-pass  channel  patterns  shown  in  Fig.  3(a)  gen¬ 
erally  result  in  a  smaller  path-length  difference  than  that  of  the 
2-pass  or  3-pass  patterns,  and  thus,  are  not  suitable  for  maximizing 
the  under-rib  convection.  Therefore,  the  4-  pass  and  5-pass  channel 
patterns  depicted  in  Fig.  3(a)  were  not  considered  for  the  generation 
and  geometrical  characterization  of  MPSFFs  in  this  study.  Never¬ 
theless  these  patterns  may  be  useful  when  it  is  necessary  to  control 
excessive  under-rib  convection  in  MPSFFs.  Two  spiral  channel  pat¬ 
terns,  denoted  as  the  5-  pass  (spiral)  and  7-pass  (spiral),  are  shown 
in  Fig.  3(a),  and  exhibit  a  close  geometrical  similarity  to  the  2-pass 
pattern.  In  fact,  the  spiral  patterns  generally  lead  to  very  high  under¬ 
rib  convection  in  MPSFFs,  similar  to  the  2-pass  pattern.  It  should  be 
noted  that  the  maximum  path-length  in  MPSFFs  tends  to  be  too 
large  as  the  pass  number  increases,  which  renders  basic  channel 
patterns  with  larger  pass  numbers  less  useful. 

Fig.  4  shows  the  distributions  of  the  path-length  difference,  A z, 
and  the  local-averaged  path-length,  zyl1,  in  parallel  flow-fields  con¬ 
structed  using  the  basic  1-pass,  2-pass,  3-pass,  5-pass  (spiral),  and 
7-pass  (spiral)  patterns.  The  z^11  in  Fig.  4(b)  is  calculated  by  averag¬ 
ing  the  path-lengths  of  infinitely  many  channels  in  they-direction. 
Fig.  4(a)  shows  that  the  1-pass  and  2-pass  patterns  set  the  lower 
and  the  upper  boundaries  between  which  A z  lies. 

Lower  boundary  :  Az  =  0.  (10) 

Upper  boundary  :  Az  =  2zmax  (- — — )  forz  < 

/  2z  \  Zmax  z  (11) 

Az  =  2zmJ-^-\  forz  >^. 

V^max  —  A  /  4 

The  parallel  flow-field  with  the  2-pass  pattern  has  the  largest 
path-length  difference,  A z,  but  also  has  a  zero  A z  at  the  mid-point 
of  the  field  (z  =  zmax/2).  The  3-pass  pattern  results  in  a  smaller  A z 
for  z  <  zmax/3  and  z  >  2zmax/3  compared  with  the  2-pass  pattern,  but 
a  larger  A z  for  zmax/2  <  z  <  2zmax/3.  That  is,  the  3-pass  pattern 
leads  to  a  more  uniform  distribution  of  A z  along  the  channel  than 
does  the  2-pass  pattern.  The  distribution  of  A z  in  the  parallel  flow- 
fields  with  the  5-pass  (spiral)  or  7-pass  (spiral)  patterns  seems  to 
approach  that  with  the  2-pass  pattern,  doing  so,  more  closely,  as 
the  pass  number  increases. 

Fig.  4(b)  indicates  that  the  2-pass  pattern  results  in  the  most  uni¬ 
form  distribution  of  zyjj  around  its  mean  value  of  zmax/2,  while  the 
1-pass  pattern  results  in  the  most  non-uniform  distribution.  The 
local-averaged  path-length,  zyl1,  is  observed  to  be  more  uniform 
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Fig.  4.  Geometrical  characterization  results  for  basic  parallel  flow-fields:  (a)  path- 
length  difference,  A z;  (b)  local-averaged  path-length,  z. 


(a) 


(c) 


Master  flow  field 


Division  into  passes 


3-pass  (single)  serpentine  flow  fields 


Fig.  5.  Design  procedure  for  generating  single-path  MPSFFs:  (a)  serpentine  master 
field  is  designed  first,  (b)  which  is  divided  into  several  passes  according  to  pass 
number  (npass  =  3  for  a  3-pass  pattern).  Finally,  (c)  proper  terminals  of  separate  passes 
are  connected  to  construct  a  single-path  MPSFF  with  (d)  inner-  or  (e)  edge-located 
inlet/outlet  ports. 


when  the  pattern  is  changed  from  the  3-pass  to  the  5-pass  (spi¬ 
ral)  and  7-pass  (spiral)  pattern.  By  contrast  the  differences  in  zyl1 
between  all  of  the  multi-pass  channel  patterns  are  not  noticeable 
in  Fig.  4(b). 

3.2.  Single-path  MPSFFs  for  small  cells 


A  simple,  step-by-step  design  procedure  to  generate  single-path 
multi-pass  serpentine  flow-fields  (MPSFFs)  is  presented  in  Fig.  5.  A 
master  field,  through  which  a  selected  channel  pattern  is  directed,  is 
designed  first,  as  shown  in  Fig.  5(a).  The  unit  tile  length  of  the  master 
field,  /tile,  is  predetermined  by  multiplying  the  pass  number,  npass, 
of  the  selected  channel  pattern  by  lunit  (/tiie  =  nPass  x  lunit).  Then,  as 
shown  in  Fig.  5(b),  the  master  field  is  divided  into  several  passes 
according  to  npaSs  (npass  =  3  in  the  Figure).  The  flow-field  obtained 
at  this  stage  is  similar  to  the  conventional  parallel  serpentine  flow- 
fields.  As  a  final  step,  proper  terminals  of  those  separated  passes  are 
connected  to  form  a  single  long  MPSFF,  as  shown  in  Fig.  5(c).  Note 
that  two  different  MPSFFs  can  be  generated,  as  shown  in  Fig.  5(d) 
and  (e),  wherein  the  positions  of  the  inlet/outlet  ports  are  different 
(edge  or  inner). 

Six  flow-field  designs  generated  for  a  square  active  area  of 
9  cm2  (3  cm  x  3  cm)  are  shown  in  Fig.  6,  wherein  the  channel  cells 
are  shaded  according  to  their  path-length  values.  A  total  of  225 
(15  x  15)  channel  cells  with  Iunit  =  2mm  were  needed  to  cover  the 
9  cm2  cell  area.  Fig.  6  shows  a  conventional  serpentine  flow-field 
(small  a),  a  conventional  spiral  flow-field  (small  f),  and  four  MPSFFs 
(small  b,  c,  d,  and  e)  that  have  been  generated  according  to  the  pro¬ 
posed  design  method.  The  master  fields  and  the  selected  channel 
patterns  are  also  provided  in  Fig.  6. 


Inlet  Outlet 

0.0  0.5  1.0 


Fig.  6.  Single  serpentine  flow-fields  generated  for  small  cells  with  a  9  cm2  active  area 
(3x3  cm):  (a)  conventional  serpentine  flow-field,  (b)  MPSFF  with  a  2-pass  pattern, 
(c)  and  (d)  MPSFFs  with  a  3-pass  pattern,  (e)  MPSFF  with  5-pass  (spiral)  pattern,  and 
(f)  conventional  spiral  flow-field. 
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The  MPSFF  shown  in  Fig.  6(b)  is  similar  to  the  flow-field  design 
proposed  for  coolant  flow-fields  in  prior  studies  [19,20]  to  provide 
better  cooling  of  PEMFCs.  Similarly,  the  MPSFF  shown  in  Fig.  6(c)  is 
almost  the  same  as  the  convection  enhanced  serpentine  flow-field 
(CESFF)  proposed  by  Xu  and  Zhao  [18]  to  improve  the  performance 
and  stability  of  DMFCs.  The  MPSFFs  (small  b,  c,  d,  and  e)  and  the  spi¬ 
ral  flow- field  (small  f)  in  Fig.  6  share  a  geometrical  characteristic, 
i.e.,  the  proximal  location  or  the  interlaced  pattern  of  channels  with 
dramatically  different  path-lengths.  On  the  other  hand,  neighbor¬ 
ing  channels  generally  have  similar  path-lengths  in  the  serpentine 
flow-field  (small  a),  as  indicated  in  Fig.  6(a). 

Fig.  7  shows  the  path-length  difference,  A z,  along  the  chan¬ 
nel  for  the  flow-fields  shown  in  Fig.  6.  The  maximum  path-length, 
zmax,  is  45  cm  (15  x  3  cm)  for  the  flow-fields.  The  shaded  triangular 
area  in  Fig.  7  indicates  the  regions  between  the  lower  and  upper 
boundaries  of  A z  identified  in  Fig.  4(a).  A  relatively  small  A z  for 
the  serpentine  flow-field  (small  a)  can  be  observed  in  Fig.  7,  which 
indicates  a  low  under-rib  convection  in  that  flow-field.  Conversely, 
the  MPSFFs  (small  b,  c,  d,  and  e)  are  observed  to  have  much  larger 
values  of  A z,  indicating  higher  under-rib  convection  in  these  flow- 
fields.  The  spiral  flow- field  (small  f)  has  the  largest  A z,  and  thus, 
the  under-rib  convection  in  this  flow-field  is  expected  to  be  the 
largest. 

The  path-length  difference,  A z,  in  the  MPSFF  with  the  2-pass 
pattern  (small  b)  is  generally  larger  than  that  in  the  MPSFFs  with 
the  3-pass  pattern  (small  c  and  d)  in  Fig.  7,  but  becomes  smaller 
for  zmax/3  <  z  <  2zmax/3.  Thus,  the  under-rib  convection  in  the 
MPSFFs  with  the  3-pass  pattern  is  believed  to  be  more  uniform 
compared  with  the  MPSFF  with  the  2-pass  pattern.  The  MPSFF  with 


Fig.  7.  Variation  of  path-length  difference,  A z,  along  channel  in  single  serpentine 
flow-fields:  (a)  for  flow-fields  of  small  a,  b,  and  c,  and  (b)  for  flow-fields  of  small  d, 
e,  and  f. 
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Fig.  8.  Spatial  distribution  of  local-averaged  path-length,  z,  in  single  serpentine 
flow-fields:  (a)  a  conventional  serpentine  flow-field,  (b)  MPSFF  with  2-pass  pattern, 
(c)  and  (d)  MPSFFs  with  3-pass  pattern,  (e)  MPSFF  with  5-pass  (spiral)  pattern,  and 
(f)  conventional  spiral  flow- field. 


the  5-pass  pattern  (small  e)  and  the  spiral  flow-field  (small  f)  have 
very  high  values  of  A z,  but  also  relatively  small  values  of  A z  near 
the  mid-point  of  z  =  zmax/2.  This  trend  is  quite  similar  to  the  dis¬ 
tribution  of  A z  for  the  parallel  flow-field  with  the  2-pass  pattern 
shown  in  Fig.  4(a).  The  most  non-uniform  under-rib  convection  is 
believed  to  occur  in  the  spiral  flow-field  (small  f). 

The  distributions  of  the  local-averaged  path-length,  z,  for  the 
single-path  MPSFFs  are  presented  in  Fig.  8.  The  serpentine  flow- 
field  (small  a)  shows  the  most  non-uniform  distribution  ofz,  while 
the  spiral  flow-field  (small  f)  shows  the  most  uniform  distribution. 
Note  that  the  distribution  of  z  is  generally  related  to  the  reactant 
concentration  distribution  in  the  PEMFCs,  and  also  to  the  temper¬ 
ature  distribution  in  the  cooling  plates  of  the  PEMFCs.  Fig.  8  shows 
that  regions  of  high  and  low  z  are  closely-interlaced  in  the  MPSFFs 
(small  b,  c,  d,  and  e),  thereby  leading  to  a  more  uniform  distribution 
ofz. 

The  geometrical  characterization  indices  for  the  single-path 
MPSFFs  shown  in  Figs.  6-8  are  summarized  in  Table  2.  The  aver¬ 
age  pass-length  difference,  ( Az)avg>i /3 ,  for  0  <  z  <  zmax/3  is  only  about 
5.4  cm  for  the  conventional  serpentine  flow-field  (small  a),  while  it 
is  in  the  range  of  17.0-32.9  cm  for  the  MPSFFs  (small  b,  c,  d,  and  e). 
The  spiral  flow-field  (small  f)  has  the  largest  (Az)avg  l/3  of  approx¬ 
imately  44.2  cm.  The  under-rib  convection  intensities,  q^lqtot,  are 
also  provided  in  Table  2,  which  are  about  4.6%  for  the  serpentine 
flow-field,  14.5-28.1%  for  the  MPSFFs,  and  37.7%  for  the  spiral  flow- 
field.  These  values  should  only  be  considered  as  estimates  because 
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Table  2 

Geometrical  characterization  indices  of  single  serpentine  flow-fields  generated  for  a  small  active  area  (3x3  cm). 


of  the  many  simplifying  assumptions  employed  in  the  derivation 
of  Eq.  (9).  Nevertheless,  qrib/qtot  in  Table  2  clearly  indicates  that 
the  MPSFFs  can  greatly  enhance  under-rib  convection  over  conven¬ 
tional  serpentine  flow-fields. 

The  standard  deviation,  o^,  of  the  local-averaged  path-length,  z 
(whose  (z)avg  is  22.5  cm),  is  also  provided  in  Table  2.  This  value  is 
related  to  the  spatial  uniformity  of  the  operational  conditions,  as 
explained  above.  The  serpentine  flow-field  (small  a)  has  the  largest 
deviation  of  about  12.7  cm,  while  the  spiral  flow-field  (small  f)  has 
the  smallest  of  approximately  3.6  cm.  The  MPSFFs  (small  b,  c,  d, 
and  e)  have  o ^  in  the  range  of  5.3-8.4cm,  which  indicates  better 
uniformity  in  operating  conditions  in  the  MPSFFs  compared  with  a 
conventional  serpentine  flow-field. 

3.3.  Parallel-path  MPSFFs  for  large  cells 

A  step-by-step  procedure  to  design  multi-pass  serpentine  flow- 
fields  (MPSFFs)  having  parallel  paths  is  presented  in  Fig.  9.  The 
similarity  between  Fig.  5  for  the  design  of  single-path  MPSFFs 
and  Fig.  9  for  parallel-path  MPSFFs  should  be  noted.  A  master 
field  shown  in  Fig.  9(a)  is  first  divided  into  several  paths  accord¬ 
ing  to  the  number  of  parallel  paths,  npath  (npath  =  3  in  the  Figure), 
then  further  into  several  passes  according  to  npass  (also  npas s  =  3  in 
the  Figure).  The  intermediate  flow-field  shown  in  Fig.  9(c)  is  the 
same  as  the  conventional  serpentine  flow-field  with  nine  paral¬ 
lel  paths.  Finally,  proper  terminals  of  those  separated  passes  are 
connected  to  form  a  MPSFF  having  parallel  paths,  as  shown  in 


Master  flow  field  Division  into  paths 


Division  into  passes  &  3-pass  (parallel)  serpentine 

Connection  of  pass  ends  flow  field 


Fig.  9(d).  The  unit  tile  length  of  the  master  field,  Ztile,  is  determined 

as  /tile  =  ft  path  x  ftpass  x  /unit  • 

Fig.  10  shows  six  flow-field  designs  for  a  square  cell  area  of 
81  cm2  (9  cm  x  9  cm),  with  three  parallel  paths  (large  a,  b,  and  c) 


Fig.  9.  Design  procedure  for  generating  parallel-path  MPSFFs:  (a)  serpentine  master 
field  is  designed  first,  which  is  divided  into  several  passes  according  to  (b)  path 
number  (npath  =  3  for  three  parallel  paths),  and  then  according  to  (c)  pass  number 
(nPass=3  for  3-pass  patterns).  Finally,  (d)  proper  terminals  of  separate  passes  are 
connected  to  construct  a  parallel-path  MPSFF. 


Fig.  10.  Parallel  serpentine  flow-fields  generated  for  large  cells  with  81  cm2  active 
area  (9  x  9  cm):  (a)  conventional  3-path  serpentine  flow-field,  (b)  3-path  MPSFF  with 
3-pass  patterns,  (c)  3-path  MPSFF  with  5-pass  patterns,  (d)  conventional  5-path 
serpentine  flow-field,  (e)  5-path  MPSFF  with  3-pass  patterns,  (f)  5-path  MPSFF  with 
9-pass  spiral  patterns. 
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or  five  parallel  paths  (large  d,  e,  and  f).  A  total  of  2025  (45  x  45) 
channel  cells  with  !unit  =  2mm  were  required  to  construct  these 
flow-fields.  The  maximum  path-length,  zmax,  was  135  cm  for  the 
3-path  parallel  flow-fields  (large  a,  b,  and  c)  and  81  cm  for  the  5- 
path  parallel  flow- fields  (large  d,  e,  and  f).  The  total  path-length, 
summing  all  of  the  parallel  paths,  equals  405  cm  (45  x  9  cm).  The 
flow-fields  shown  in  Fig.  10  can  be  categorized  as  conventional  par¬ 
allel  serpentine  flow-fields  (large  a  and  d),  and  parallel-path  MPSFFs 
constructed  based  on  a  3-pass  channel  pattern  (large  b  and  e),  a 
5-pass  spiral  pattern  (large  c)  and  a  9-pass  spiral  pattern  (large  f). 
Highly-interlaced  channel  patterns  in  the  parallel-path  MPSFFs  can 
clearly  be  observed  in  Fig.  10. 

The  distribution  of  the  path-length  difference,  A z,  along  the 
channel  is  illustrated  in  Fig.  11.  The  serpentine  flow-fields  (large 
a  and  d)  have  very  small  values  of  A z.  The  reason  for  this  trend 
is  because  most  channels  are  in  contact  with  channels  of  similar 
path-lengths  due  to  the  parallel  arrangement  (see  Figs.  10(a)  and 
(d)).  The  MPSFFs  (large  b,  c,  e,  and  f)  have  much  larger  values  of  A z 
compared  with  the  serpentine  flow-fields,  which  suggests  that  the 
under-rib  convection  in  these  MPSFFs  is  greatly  enhanced.  Fig.  11 
also  shows  that  the  3-pass  pattern  leads  to  a  more  uniform  distri¬ 
bution  of  A z  than  does  the  5-pass  spiral  pattern  or  the  9-pass  spiral 
pattern.  A  uniform  distribution  of  A z  is  important  for  obtaining  a 
uniform  under-rib  convection  intensity  throughout  the  flow-field. 
The  similarity  between  the  spiral  channel  patterns  and  the  2-pass 
pattern  is,  again,  demonstrated  in  Fig.  11,  i.e.,  A z  of  the  MPSFF  with 
the  9-pass  spiral  pattern  closely  follows  the  upper  boundary  that 
has  been  determined  from  the  parallel  flow-field  with  the  2-pass 
pattern  in  Fig.  4(a). 


Fig.  11.  Variation  of  path-length  difference,  A z,  along  channel  in  parallel  serpentine 
flow-fields:  (a)  for  3-path  flow-fields  of  large  a,  b,  and  c,  and  (b)  for  5-path  flow-fields 
of  large  d,  e,  and  f. 


Fig.  12.  Spatial  distribution  of  local-averaged  path-length,  z,  in  parallel  serpentine 
flow-fields:  (a)  conventional  3-path  serpentine  flow-field,  (b)  3-path  MPSFF  with 
3-pass  patterns,  (c)  3-path  MPSFF  with  5-pass  patterns,  (d)  conventional  5-path 
serpentine  flow-field,  (e)  5-path  MPSFF  with  3-pass  patterns,  (f)  5-path  MPSFF  with 
9-pass  spiral  patterns. 


Fig.  12  shows  that  the  local-averaged  path-length,  z,  is  more  uni¬ 
formly  distributed  in  the  MPSFFs  (large  b,  c,  e,  and  f)  than  in  the 
serpentine  flow-fields  (large  a  and  d),  indicating  that  the  operation 
conditions  are  expected  to  be  more  uniform  in  the  MPSFFs.  This  also 
points  out  that  these  parallel-path  MPSFFs  can  work  well  as  coolant 
flow-fields  for  PEMFCs  or  DMFCs.  Fig.  12  also  shows  that  increasing 
the  parallel  path  number,  npath,  enhances  the  uniformity  of  the  z 
distribution  in  the  cell  area. 

The  geometrical  characterization  indices  for  the  parallel-path 
MPSFFs  shown  in  Figs.  10-12  are  summarized  in  Table  3.  The 
MPSFFs  (large  b,  c,  e,  and  f)  are  observed  to  have  significantly 
larger  average  path-length  differences,  (Az)avg)i/3,  compared  with 
the  serpentine  flow-fields  (large  a  and  d).  The  under-rib  convection 
intensity,  qrib,i/3/<7tot,  is  only  about  5.4%  or  2.9%  for  the  serpentine 
flow-fields  (large  a  and  d),  while  it  is  278.3%  or  109.4%  for  the 
MPSFFs  with  the  3-pass  pattern  (large  b  and  e).  Although  a  value 
f°r  9rib,i/3/<7tot  greater  than  100%  definitely  violates  the  assumption 
of  a  negligible  under-rib  convection  for  Eq.  (9),  such  large  values 
of  Qrib.i/3/^tot  clearly  point  to  significantly  higher  under-rib  convec¬ 
tions  in  those  parallel-path  MPSFFs. 

Table  3  indicates  that  an  increase  in  the  path  number,  npath,  from 
three  (large  a,  b,  and  c)  to  five  (large  d,  e,  and  f)  generally  decreases 
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Table  3 

Geometrical  characterization  indices  of  parallel  serpentine  flow-fields  generated  for  a  large  active  area  (9x9  cm). 


the  under-rib  convection  intensity.  This  is  explained  by  the  linear 
proportionality  of  qrib,i/3/9tot  to  zmax  given  in  Eq.  (9).  The  maxi¬ 
mum  path-length,  zmax,  is  inversely  proportional  to  npath,  and  so 
is  the  under-rib  convection  intensity,  grib,i/3/9tot.  Comparison  of  the 
MPSFF  with  the  3-pass  pattern  (large  b)  and  with  the  5-pass  spi¬ 
ral  pattern  (large  c)  reveals  that  spiral  patterns  generally  result  in 
higher  under-rib  convection  intensities,  which  is  consistent  with 
the  result  of  Fig.  4(a). 

Table  3  also  suggests  that  the  standard  deviation,  oz,  of  the  z 
distribution  in  the  active  area  is  significantly  reduced  from  about 
38.8  or  23.2  cm  in  the  serpentine  flow-fields  (large  a  and  d)  to  about 
13.2  or  6.1  cm  in  the  MPSFFs  with  the  3-pass  pattern  (large  b  and 
e).  The  reduction  of  o-z  by  a  factor  of  approximately  1/3  indicates 
that  parallel-path  MPSFFs  are  favourable  for  the  uniformity  of  the 
local-averaged  path-length,  z,  and  thus  also  for  the  uniformity  of 
operating  conditions.  The  inverse  proportionality  of  oz  to  the  num¬ 
ber  of  parallel  paths,  npath,  can  also  be  observed  in  Table  3. 

4.  Discussion 

This  study  focuses  on  increasing  the  path-length  difference,  A z, 
in  serpentine  flow-fields  with  the  hypothesis  that  an  enhancement 
of  under-rib  convection  between  neighbouring  channels  improves 
the  performance  of  PEMFCs.  The  results  of  Xu  and  Zhao  [18]  and  Qi 
and  Kaufman  [22]  support  the  favourable  effects  of  under-rib  con¬ 
vection  enhancement  on  PEMFC  performance.  On  the  other  hand 
the  negative  effects  of  enhanced  under-rib  convection  have  also 
been  addressed  by  Rock  [27],  who  proposed  mirrored  serpentine 
flow-field  designs  for  fuel  cells  and  suggested  that  the  performance 
of  PEMFCs  could  be  improved  by  reducing  the  cross-leakage  flow 
or  bypass  flow  in  the  flow-fields.  The  terms  cross-leakage  and  by¬ 
pass  flow  denote  the  same  process  as  under-rib  convection.  These 
seemingly  controversial  results  concerning  under-rib  convection  in 
serpentine  flow-fields  may  be  explained  by  the  difference  in  the  cell 
sizes  considered  in  those  studies  [18,22,27]. 

The  path-length  difference,  Az,  of  conventional  serpentine 
flow-fields  is  directly  determined  by  the  cell  size.  For  example,  a 
serpentine  flow-field  for  an  Lem  x  Lem  active  area  has  an  average 
path-length  difference,  ( Az)avg,  of  21.  Thus,  ( Az)avg  in  the  serpen¬ 
tine  flow-field  is  only  about  6  cm  for  small  cells  having  a  3  cm  x  3  cm 
active  area.  Recalling  that  the  under-rib  convection  intensity  is  pro¬ 
portional  to  the  path-length  difference,  it  is  well  understood  that 
small  cells,  such  as  the  one  considered  by  Xu  and  Zhao  [18],  can  per¬ 
form  better  when  the  under- rib  convection  is  enhanced;  however, 
large  cells,  such  as  the  20  cm  x  24  cm  cell  considered  in  Shimpalee 
et  al.  [28]  (based  on  the  design  of  Rock  [27]),  have  a  large  ( Az)avg 
of  about  40  cm,  and  thus,  already  have  a  high  under-rib  convection 
intensity.  As  pointed  out  in  Rock  [27],  excessive  under-rib  convec¬ 
tion  in  serpentine  flow-fields  can  result  in  premature  exhaust  of 
reactant  gases  due  to  the  flow  short-circuiting.  This  bypass  flow 
through  the  under-rib  region  also  reduces  the  flow  speed  in  the 
channel,  which  adversely  impacts  reactant  distribution  and  liquid 
water  transport.  In  this  case,  higher  performance  is  achieved  by 
regulating  the  excessive  under-rib  convection  in  PEMFCs. 

The  MPSFFs  proposed  in  this  study  provide  a  way  to  maximize 
the  path-length  difference,  and  thus,  the  under-rib  convection  in  a 


given  cell  area.  Based  on  the  above  discussion,  MPSFFs  are  expected 
to  be  more  advantageous  for  small  fuel  cells  primarily  employed 
for  portable  applications.  The  enhanced  under-rib  convection  in 
MPSFFs  decreases  the  pressure  drop  across  PEMFCs,  which  also 
contributes  to  the  performance  by  reducing  the  power  consump¬ 
tion  of  air  blowers  in  portable  systems.  More  uniform  distributions 
of  reactant  concentration,  temperature  and  liquid  water  saturation 
are  another  advantage  of  MPSFFs,  which  enable  better  control  of 
PEMFC  operation,  as  well  as  the  use  of  dry  reactant  gases  with  inter¬ 
nal  humidification.  The  MPSFFs  are  also  expected  to  be  useful  to 
achieve  better  uniformity  in  coolant  flow-fields  for  PEMFCs. 

This  study  geometrically  characterized  single-  and  parallel-path 
MPSFFs  based  on  the  path-length,  z,  and  estimated  the  under-rib 
convection  intensities  in  the  flow-fields.  Nevertheless,  the  exact 
effects  of  the  enhanced  under-rib  convection  on  the  performance 
of  PEMFCs  still  need  to  be  clarified  through  further  numerical 
and  experimental  investigations.  For  this  purpose,  the  authors  are 
currently  working  on  numerical  studies  to  determine  the  flow  dis¬ 
tribution,  pressure  drop,  and  cooling  characteristics  of  the  proposed 
MPSFFs. 

5.  Summary 

A  simple,  step-by-step  design  method  to  generate  single-  and 
parallel-path  multi-pass  serpentine  flow-fields  (MPSFFs)  for  PEMFC 
applications  is  presented  in  this  study.  MPSFFs  comprised  of  multi¬ 
pass  flow  channels  arranged  in  serpentine  configurations  allow 
channels  with  widely  different  path-lengths  to  be  positioned 
near  one  another.  The  resulting  maximum  path-length  differ¬ 
ence  between  neighbouring  flow-channels  is  expected  to  enhance 
under-rib  convection  and  transport,  thereby  improving  the  per¬ 
formance  of  PEMFCs.  In  addition,  the  highly-interlaced  channel 
patterns  in  MPSFFs  are  expected  to  improve  the  uniformity  of  local 
conditions,  such  as  reactant  and  product  concentrations,  temper¬ 
ature  and  liquid  water  saturation  in  the  active  cell  area.  Several 
single-  and  parallel-path  MPSFFs  generated  by  the  proposed  design 
method  have  also  been  presented,  along  with  the  results  of  their 
geometrical  characterization.  The  implication  of  under-rib  convec¬ 
tion  enhancement  on  the  performance  of  PEMFCs  with  different 
cell  sizes  has  also  been  discussed. 
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